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ABSTRACT 

We apply the supernova(SN) extinction curves to reproduce the observed properties 
of SST J1604+4304 which is a young infrared (IR) galaxy at z ~ 1. The SN extinction 
curves used in this work were obtained from models of unmixed ejecta of type II 
supernovae(SNe II) for the Salpeter initial mass function (IMF) with a mass range 
from 8 to 30 Mq or 8 to 40 Mq. 

The effect of dust distributions on the attenuation of starlight is investigated 
by performing the % 2 fitting method against various dust distributions. These are the 
commonly used uniform dust screen, the clumpy dust screen, and the internal dust ge- 
ometry. We add to these geometries three scattering properties, namely, no-scattering, 
isotropic scattering, and forward-only scattering. Judging from the x 2 values, we find 
that the uniform screen models with any scattering property provide good approxi- 
mations to the real dust geometry. Internal dust is inefficient to attenuate starlight 
and thus cannot be the dominant source of the extinction. 

We show that the SN extinction curves reproduce the data of SST J1604+4304 
comparable to or better than the Calzetti extinction curve. The Milky Way extinction 
curve is not in satisfactory agreement with the data unless several dusty clumps are in 
the line of sight. This trend may be explained by the abundance of SN-origin dust in 
these galaxies; SN dust is the most abundant in the young IR galaxy at z ~ 1, abundant 
in local starbursts, and less abundant in the Galaxy. If dust in SST J1604+4304 is 
dominated by SN dust, the dust production rate is ~ 0.1 Mq/SN. 

Key words: dust, extinction — galaxies: evolution — galaxies: ISM — galaxies: 
starburst — infrared: galaxies — supernovae: general 



1 INTRODUCTION 

Dust is a minor constituent in the Universe, but it plays 
crucial roles for evolution of various objects; for example, 
it cools gas clouds by radiating far-infrared (far-IR) emis- 
sion, and acts as a catalyzer to form molecules. Correc- 
tions for dust extinction have significant cosmological con- 
sequences for determining the Hubble constant, the cosmo- 
logical model, and the epoch of galaxy formation. 

It is generally accepted that the atmospheres of asymp- 
totic giant branch(AGB) stars are one of the main sources 
of interstellar dust in the Milky Way (MW). However, AGB 
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stars are too old to account fo r the presence of dust at z > 6 
(but see Valia nte et al.l [2009D . and supernovae (SNe) have 
been recognized as a candidate of the major source of dust 
in the early stage of galaxy evolution. 

Stars with masses 8-40 Mq end up their life 
as type II supe rnovae(SNe Ilfl in 40 Myr since 
their formation (|Heger et al.l l200&t ). The formation of 
dust in SN ejecta was first witnessed in SN 1987A, 



1 In this paper, we simply use the term SNe II for core-collapsed 
SNe whose progenitors are massive stars. We also call the dust 
produced by SNe II "SN dust" , and the extinction curve synthe- 
sized for the SN dust "SN extinction curve" . 
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motivating Kozasa. Hasegawa fc Nomotol l|l989l ) and 
iKozasa, Hasegawa fc Nomotol ( 199ll ) to predict the con- 
densed dust composition and size. Systematic studies 
of dust formation in SN ejecta have been proceeded to 
explore the role of dust in the formation and evolution of 
stars and galaxies in the e arly universe dTodini fc Ferrara 
200 ll : iNozawa et all 120031 ; | Schneider. Ferrara. fc Savaterral 



2004]): the composition of synthesized dust reflects the 
elemental composition in the formation site which is 
affected by the type and degree of mixing and the degree 
of formation efficiency of CO molecule, and the size of 
dust is controlled by the time evolution of g as density and 
temperature depending on the type of SNe |Nozawa et al.l 
l2008llNozawa et al.ll201ol ). 

Apart from the differences in the SN models and the 
underlying assumptions employed for dust formation cal- 
culations, the predicted dust mass formed in the ejecta 
reaches 0.1 to IMq per SN II, which is sufficient enough 
to explain a vast amount of dust of > lO 8 A/0 discov- 
ered in host g alaxies of high-redsh i fted q uasars (QSOs) 
at z > 5feg. iMorgan fc Edmunds! 120031 ; IMaiolino et all 
120061 ; iDwek.. Galliano fc Jones! 120071) . However, the mass 
and composition of dust formed in the ejecta have been still 
controversial; the observations of nearby SNe have claimed 
dus t mass formed in t he ejecta to be less than 1O _3 M0 
fsee lKozasa et al.ll2009l for a review), while the infrared and 
submillimeter observations of Galactic supernova remnants 
(SNRs) h ave suggested the presence of dust reaching 0.02 



to l Mfl jRho et al.ll20qsl; iDunne et al.ll2003l; iDunne et al 



20091 ; iGomez et al.ll2009l ; ITemim et al.ll2010UNozawa eta! 
2010). Although no information has been available on the 
compos ition of dus t formed in the ejecta except for SN 
2004et (Kotak 20091) , the spectral mapping observation cov- 
ering the entire of Cas A SNR with the Spitzer Space Tele- 
scope have revealed the presence of a variety of dust species 
whose composition is closely associated with the gas emis- 
sion lin e characteristic o f nucleosynthetic layer in the pro- 
genitor jRho et al.ll2003 ). 

The amount, composition and size of dust supplied from 
SNe to the interstellar medium (ISM) are not always the 
same as those condensed in the ejecta due to the processing 
in a hot plasma produced by the reverse and forward shocks 
during the journey to the ISM jBianchi fc S chneider! 120071 ; 
INozawa et alj 20071) as well as in the ho t / warm region of ISM 
(I Jones et all [l994l; Ijones et all [l996l; INozawa et al l 120061 ; 
iHirashita et all |2010| ). Nevertheless IMaiolino et all (|2004l ) 
found that the extinction curve of SDSS J1048+4647, a 
broad absorption line (BAL) QSO at z~6.2 is different 
from those of nearby galaxies, and is in excellent agreement 
with the extinction curv e obtained from SN dust models by 
iTodini fc Ferraral ([2001) considering the uniformly-mixed 
elemental composition in the ejecta. Synthes izing dust ex- 
tinction curves for various SN dust models bv lNozawa et al] 
l|2003l ). IHirashita et all (|2005l ) have shown that the extinc- 
tion curve observed in this BAL QSO can be reasonably 
reproduced by the dust produced in SNe II without any 
mixing in the ejecta. Also the extinction curve derived from 
the analysis of the afterglow of 7~ray burst GRB 050904 at 
2 = 6.3 have suggested that SNe are the domin ant source 
of dust in t he ho st galaxy (|Stratta et al.l 12007). although 
IZafar et al.l l|2010l ) have claimed no evidence of dust extinc- 
tion in the afterglow. However, the observed extinction curve 



in the local environments may provide less constraint on the 
SN dust model as well as the role of SNe on the evolution 
of galaxies. 

Young galaxies, in which no AGB stars evolved off the 
main sequence, are obvious targets to study the properties 
of SN dust, for example, extinction curves that determine 
the relation between changes in color and total absorption. 
However, it is generally difficult to infer the amount of the 
dust obscuration, because different geometries of the dust 
spacial distribution result in different observed colors for 
the same amount of dust and the same i ntrinsic spectrum 
of the underlying stellar population(e.g., iGiavaliscq |2002| ) . 
In fact, some studies of HI recombination lines in starburst 
and dusty Seyfert galaxies suggest that the extinction ob- 
served in the near-IR is ge nerally greater than that observed 
by using the optical lin es (|Kawara. Nishida fc Phillipslll989l ; 
iPuxlev fc Brandll 19941). 

ICalzetti. Kinney, fc Storchi-Bergmannl (1 19941 ) has 
shown that the extinction laws in the Milky Way and Large 
Magellanic Cloud cannot account for the ultraviolet (UV) to 
optical spectra of UV-bright starburst galaxies in the local 
Universe in a satisfactory way. They empirically derived 
the extinction law for these galaxies. Their extinction law 
is characterized by an overall slope which is more gray 
or flatter than the MW extinction law's slope and by the 
absence of the 2175 A dust feature. In starburst galaxies, 
numerous SN II explosions occur and SN-condensed dust 
would be injected into the ISM. Thus, newly formed SN 
dust would produce a great impact on the extinction law in 
local starburts. 

In this paper, we apply synthetic exti nction curves on 
the u nmixed ejecta models of SNe II by [Hirashit a et al.l 
(2005) to the broad-band SED of a young ultra luminous 
IR galaxy (ULIRG) S ST J1604+4304 at z ~ 1. Compar- 
ing with the MW and ICalzettil (|200ll ) extinction laws, we 
show that the SN extinction curves fit better to the obser- 
vations. We adopt Ho — 70 km s _1 Mpc -1 , Q. m — 0.3, and 
SIa = 1 — throughout this paper. 



2 DATA 

Synthetic spectra of stellar populations reddened by dust are 
modelled to reproduce the data of SST J1604+4304. This is 
a young IR galaxy at z — 1.135 with the characteristic of 
ULIRGs. The multi- wavelength data from the X-ray to radio 
are compiled in lKawara et al.l (|2010l ). There are no support- 
ive evidence for AGN activities in this galaxy, although the 
data are not deep enough to rule out the presence of an 
AGN. 

Figure [T] illustrates the broad-band SED from the 
UV to far-IB^ along with representative SED templates 
of starburst-dominated ULIRGs, namely, Arp220(local 
ULIRG) and MIPS J142824.0+352619 (hyper-luminous IR 
galaxy at z = 1.3). The broad-band SED and stellar ab- 
sorption lines, such as Call and HI lines, as well as the 
[OII]A3727 line suggest that young stellar populations dom- 
inate the optical luminosity in this galaxy. The Call index, 



2 The data are taken from Table 1 in iKawara et al" and 
the correction for the gravitational amplification is applied. 
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Figure 1. UV to far-IR SED for SST J1604+4304 com- 
pared with those of star burst-dominated ULIRGs, Arp 
220(solid local ULIRG), MIPS J142824.0+352619(das/i z = 
1.3 hyperluminous IR galaxy). Arp 220 data come from 
ICarico et al.l |l990l). iR igopo ulou, Lawrence. fc Rowan- Robinsonl 



ll!996ll. and iKlaas et al.l Jl997tl. and MIPS J 142824.0+352619 
from lDesai et all ^OOrf and lBorvs et alj d2006l) . 



defined by iKawara et all (j2010l ). is D(CaII) = 0.03 ± 0.11 
and the equivalent [Oil] width is EWQOII]) = 31.1 ± 4.3 
A in the rest-frame which corresponds to the observed [Oil] 
line flrafl of 1.9 ± 0. 3 x 1 0~ 17 ergs cm -2 s _1 . 

IKawara et al.l (|2010T ) estimated the reddening to be 
E(B — V) ~ 0.8 when the Calzetti extinction law is ap- 
plied. The 8-1000um bolometric luminosity inferred from 
the three Spitzer MIPS bands at the 24, 70, and 160um 
is Li r = 1.78 ± 0.63 x 10 12 Lq, corresponding to the flux 
F ir = 1.00 ± 0.35 x 10~ 12 ergs cm" 2 s _1 . These values are 
intrinsic to the galaxy after correcting for a gravitational 
amplification of 1.17 by the foreground galaxy cluster CL 
1604+4304 at z = 0.9. The bolometric luminosity agrees 
with the extinction-corrected, stellar luminosity, suggesting 
SST J1604+4304 can be characterized by models of pure 
star formation with no AGN emission. The dust mass in- 
ferred from the far-IR data is ~ 2 x 10 s Mq, sufficient to 
form a shell surrounding the galaxy and produce the ob- 
served extinction. 

The SED data at 13 photometric bands from B to 
5.8um are used for the following analysis; the longest 
Spitzer IRAC 8.0um band is not used, because the dust 
emission may be greater than the starlight in this band cor- 
responding to 3.7um in the rest-frame. 



3 EXTINCTION 

3.1 Geometry of dust distribution 

The attenuation produced by dust in front of an extended 
source does not only depend on the properties of dust grains, 
such as cross sections, albedo, and amount, but also de- 
pends on the geometry of spacial distribution of dust. It is 
known that the apparent optical depth is lower than the 



3 The [Oil] line flux published in IKawara et al.l <201Cft is incor- 
rect. 



optical depth averaged over the source area when starlight 
is attenuated by a thick, clumpy dust layer or dust is 
internal to the extended source. We consider the follow- 
ing geometries of d ust di s tribution. For full deta i ls, th e 
rea der may refer to ICodel ll 19731) iNatta fc Panagial l|l984h . 
and ICalzetti. Kinney, fc Storchi-Bergmannl 1 19941 ). 

Internal dust - In the case that absorbers (dust grains) 
and emitters(stars) are uniformly distributed, the apparent 
extinction is given by: 
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= -ln[(l - e" 
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where r ap is the apparent optical depth of the inter- 
nal dust, T e ff the effective absorption optical depth, r ex t 
the sum of the optical depth for absorption and scat- 
tering, 7 the albedo, /° is the unobscured, intrinsic flux 
density , and f v is the reddened, observed flux density. 
g =< cos{9) > is the scattering asymmetry factor, where 
g — 1 for the forward-only scattering and g = for 
the isotropic scattering. Sometimes r e / / = (1 — 'y)°' 5 T ex t 
is us e d for the isotropic scattering (i.e., Natta fc Panagial 



1 1984 ICalzetti. Kinney, fc Storchi-Bergmannl 1 19941 ) ~ How- 
ever, r e ff = (1 — j) Text approximates better to the solu- 
tion by the two-stream approximation for the geometry with 
a uniform distribution of the absorbers and the emitters. For 
the full details, refer to Figure lATl in Appendix A. 

Dust screen - In this geometry, the absorbers are phys- 
ically separated from and located in a screen in front of 
the emitters. We adopt an analytic expression derived by 
iNatta fc Panagial (|l984l ) for a system of clumps distributed 
according to Poisson statistics. 
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where N is the average number of clumps along the line 
of sight and all clumps are assumed to have one and the 
same effective optical depth r e ///JV. For N — ¥ 00, r ap — > 
T e ff, corresponding to the uniform dust screen. For the dust 
screen which is physically distant from the system of stars, 
the effect of the dust is to remove photons from the line of 
sight through absorption and scattering. On the other hand, 
for the dust screen which is located close to the system of 
stars, the scattered photons can get in the line of sight. Thus, 
T e ff is given as: 



Teff = r ext 

T ^ff = (1 - j)Text 



for no scattering (7) 
for scattering with g — 1 (8) 



Teff = (1 — 7) ' r ex t for scattering with g = 



(9) 



Again, g = 1 is for the forward-only scattering, and g = 
for the isotropic scattering. 

Composite geometry - It is quite likely that dust grains 
distribute within the system of stars as internal dust and in 
the screen in front of the system. We approximate this case 
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Figure 2. The apparent optical depth r ap is shown as a func- 
tion of the effective optical depth T e f f . The solid lines represent 
dust screen geometries for N = [1,2,4,8,16,32,64,128,256, 512,oo] 
from bottom to top. The dashed line represents an internal dust 
geometry. 



Figure 3. Syntheti c SN extinct ion curve compared wi th those 
for th e Milky way l|Drainell2003r i and local starbursts ijCalzettil 
l200ll) . The curves are normalized to the extinction (absorption + 
scattering) in magnitude at 0.8 um. Note that a spike on the SN 
curve is produced by crystal Si02 at 8.7 um. 



Figure[2]shows the apparent optical depth r ap as a func- 
tion of r e ff. It is apparent that the growth of r ap turns over 
at a smaller value of r e // in the dust screen with a smaller 
value of N. As shown in this figure, the extinction by dust 
in the screen with a small N or by internal dust strongly 
depends on the geometry and tends to produce a gray ex- 
tinction, and dust in a uniform screen produces the most 
effective extinction. 



3.2 Extinction curves 

We examine two extinction curves synthesized in eject a of 
SNe I I, the empirical extinction curve derived by ICalzett ' 



2001), and the MW extinction curve given by 



Drains 



2003). Figure [3] plots these extinction curves as a function 



of wavelengths, showing that the SN extinction curves are 
more gray or flatter than the Calzetti and MW extinction 
curves from 0.1 - lum. 

The MW extinction curve exhibits a pronounced fea- 
ture at 2175 A. This feature is absent in the typical extinc- 
tion curves of the Small Magellanic Clouds (SMC). The dust 
in local starbursts and Lyman break galaxies e xhibits the 
prope rties of the SMC- like dust properties, while [Noll et al.l 
|2009l ) claimed that about 30% of the spectra of actively 
star-forming, intermediate-mass galaxies at 1 < z < 2.5 ex- 
hibit broad absorption features at ~2175 A. To check the ef- 
fect of the feature on the extinction law in SS T J1604+4303 , 
in addition to the MW extinction curve bv lDrainej (|2003l ). 
we will examine the pseudo MW extinction curve which is 
made by removing the 2175 A feature from the MW extinc- 
tion curve. 

The ICalzettil |200ll ) extinction curve was derived by 
analyzing local starbursts and blue compact dwarf galaxies 
l|Calzetti. Kinney, fc Storchi-Bergmannll 19941 ) , and is widely 
used to, for exa mple, derive pho tometric redshifts of high- 
redshift galaxies. ICalzettil (|200lf ) gives the extinction curve 
in a narrow range from 0.12 to 2.2 um. To analyze multi- 
wavelength data of high-redshift objects from the UV to IR 
such as given in Figure [T] the extinction curve should be 



ex tended to longw ard of 2.2 urn. The functional form given 
by ICalzettil l|200ll ) cannot be used to extend the extinction 
curve, because the extinction value sharply drop beyond 2.2 
pi. The Calz etti extinction curve is almost identical to the 
iDraind (|2003h 's MW curve from 0.6 to 1.6 um. We thus ex- 
tend the Calzetti curve to longer wavelengths by adopting 
the MW extinc tion curve for w avelengths 1.6 um. In ad- 
dition, because ICalzettil l|200ll) does not pro vide the albedo, 
we adopt the albedo given in iDraine] l|2003l ) with removing 
the dip in the albedo which is attributed to the familiar 2175 
A feature. This curve is called the Calzetti extinction curve 
in this paper. It is quite likely that the Calzetti extinction 
curve reflects the dust composition from the two origins, 
namely, dust condensed in expanding atmospheres of AGB 
stars and ejecta of SNe II, because in local starburt galaxies 
low-mass stars have already evolved off the main sequence 
and reached to the dust forming stages. Because the Calzetti 
curve was empirically obtained by assuming a geometry of 
uniform foreground screens (i.e., /„ = e~ T f°), we apply this 
curve in this form only. 

The two SN extinction curves are compared with the 
MW and Calzetti extinction curves. The SN extinction 
curves were synthesiz ed in unmixed ejecta of SNe II by 
iHirashita et al.l (120051), bas e d on the dust produc tion mod- 
els bv lNozawa et all (|2003l ). iNozawa et"afl (120031 ) have in- 
vestigated the dust condensation in SNe II with various 
progenitor masses, taking into account the time evolution 
of gas temperatures in the ejecta by solving the radiative 
transfer equation with energy deposition of radioactive ele- 
ments. The two extreme cases are considered for the elemen- 
tal composition in the He-core: the unmixed case and the 
uniformly mixed case. The formation of all possible conden- 
sations is calculated simultaneously by applying a theory of 
non-steady state nucleation and grain growth under the as- 
sumption of complete formation of CO and SiO molecules. 
Various dust species condense according to the elemental 
composition at the formation site in the unmixed ejecta; 
Carbon in the inner C-rich He layer, AI2O3, Mg-silicates 
(Mg2Si04 and MgSiOs) in the oxygen-rich layer, MgO in 
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the outer O-rich layer, Si02 in the inner O-rich layer, Si 
and FeS in the Si-S layer, and FeS and Fe in the inner most 
Fe— S layer. In the mixed ejecta, only oxide grains such as 
AI2O3, Mg-silcaites, Si02 and FesCU condense because the 
elemental composition assuming complete formation of CO 
is oxygen-rich in the entire He-core. In the unmixed cases, 
C and Si dust sp ecies dominate the dust mass in the low- 
mass progenitor (jNozawa et al.ll2003l ). and Mg2Si04 is get- 
ting dominant with increasing progenitor mass, while in the 
mixed ejecta Si02 and Mg2Si04 dust dominate the dust 
mass. 

Combining the optical constants in the literature with 
the grain s pecies and size distrib ution of SN-condensed dust 
derived by Nozawa et all l|2003l ), and using the Mie theory, 
iHirashita et all (|2005h have calculated absorption and scat- 
tering cross sections of homogeneous spherical grains with 
various sizes condensed in mixed and unmixed ejecta. It 
should be pointed out here that in the unmixed ejecta C 
and Si dust grains with radius a ^ O.l^tm are significantly 
abundant in the low mass progenitors than in the high mass 
progenitors. This results in the different behavior of the ex- 
tinction curves: the extinction curve in UV region tends to 
be flatter in the unm ixed cases than in the mixed cases 
l|Hirashita et al.l I2005T ) . The details for the a dopted quan- 
tities are described in IHirashita et alj |2005l ). The extinc- 
tion curves in unmixed ejecta of SNe II reasonably repro- 
duce the extinction l aw of the z = 6.2 BAL QSO found by 
iMaiolino et al.l |2004l '). In this paper, we use the two extinc- 
tion curves of unmixed SNe II for the lSalpeterl (|l955t ) Initial 
Mass Function (IMF); one has a mass range of progenitors 
from 8 to 30 Mq, and the other from 8 to 40 Mq. These 
were obtain ed by interpolat i ng or extrapolating extinction 
curves that IHirashita et al.l (|2005t ) calculated for unmixed 
SNe II with progenitor masses of 13, 20, 25, & 30 MqQ. 



SYNTHETIC SPECTRA OF STELLAR 
POPULATIONS 



The e volutionary synthesis codes by iBruzual fe Chariot! 
(2003) are used to generate synthetic spect ra of evo l ution - 
ary stellar population models. We adopt the lSalpeterl (|l955T ) 
initial mass function (IMF) with the lower and upper mass 
cutoffs of 0.1 Mq and IOOMq, respectively. It is noted that 
the weight of massive stars is greater in the Chabrier IMF 
than in the Salpeter IMF, resulting in a mass to luminosity 
ratio roughly 1.6 times heavier in the Salpeter IMF. The star 
formation histories are exponentially declining models with 
the star formation rate (SFR) oc tj} r exp(-t/tsf r ), where t s f T 
is the e- folding timescale and t is the age, i.e., the time after 
the on set of the first star formation. Although lKawara et al.l 
(|2010h used the instantaneous burst models, such models 
are not appropriate to reproduce simultaneously the strong 
[Oil] A 3727 emission line from HII regions ionized by young, 
hot stars and the observed SED. The synthetic spectra are 
generated in a parameter range of t = 3 x 10 5 — 8 x 10 9 yr in 



4 In Hirashit a et al. the slope of extinction curve for 25 

Mq and 30 Mq progenitors was overestimated because of the 
error in the calculation processes. We have repl aced the results 
with c orrect ones. Nevertheless, the conclusions in lHirashita et ahl 
i2005l l are not altered by this change. 



a step of 0.2 dex and t s fr = 10 6 — 1.6 x 10 10 yr in a step of 
0.3 dex. For the shortest t s f r , the star formation history is 
near an instantaneous burst, while for the longest t s f r , the 
star formation rate is almost constant over the cosmic time. 
The metallicity of the synthetic spectra has a range of Z — 
0.005 - 2.5 Z Q (i.e., six points at 0.005, 0.02, 0.2, 0.4, 1.0, & 
2.5 Zq). 

The extinction models discussed in section 3 are then 
applied to the synthetic spectra generated in the above pro- 
cedure, yielding template spectra which are used for the sub- 
sequent x 2 -fitting. In the initial run, the extinction measured 
at 0.3 (xm has a range of T e:ct (0.3um) = — 10 in a step of 
0.1. In case that the initial range is too narrow to have the 
minimum \ 2 for a certain model, the range and the step are 
widen by a factor of three to search for the minimum x 2 j 
and so on. 



5 x 2 TEMPLATE-FITTING METHOD 

We search for the best set of parameters for SST J1604+4304 
by using a \ 2 template-fitting method as defined as below: 



X = Xsed + xfon] + Xcaii + X% 



(11) 



where the degree of freedom is v = 12. In what follows, we 
define the terms in equation (11). 
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where b is the parameter determined by fitting the template 
flux density Ftempiate to the observed F b s - F b s are the pho- 
tometric data measured in 13 broad band filters from B to 
5.8(ixn. The Spitzer IRAC 8.0 um is not used, because dust 
emission may be significant in this band corresponding to 
3.7 um in the rest-frame. It is noted that 5% of the flux was 
added in the quadratic form to the photometric error a p b, 
i.e., o% ED = (Tp h + (0.05F o bs) 2 ■ This is necessary to avoid 
too small and non-realistic photometric error. 

The strong [OII]A3727 line requires the presence of OB 
stars. We thus include xfon] m the X 2 template-fitting 
method, xfon] ^ defined as: 



2 

X[on] 



T [OII] 



-EW, 



[Oil] ■ 
template , 



2 / 2 
/ a [OII] 



(13) 



r \oii] 



where EW^' 1 = 31.1 ± 4.3 A. To derive EWl mplate , 
first obtain the H a luminosity by counting ionizing photons 
in the synthesized spectra. We then adopt the [OH] and H Q 
calibration performed by iKewlev. Geller fe Janser] l|2004l ). 
For the metallicity scaling of [OII]/H a , we use their theo- 
retical curves for q = 3 x 10 7 cm s _1 with an additional 
condition that [OII]/H a should not be less than 0.1. This 
additional condition is applied to avoid a negative or un- 
realistic value of [OII]/H Q outside of the calibrated range 
of metallicity. The calibrated [OII]/H a has a peak value of 
1.6 at Z = 0AZQ and dr ops to a value of 0.1 at Q14Z0 
or Z = 2.44Z0. Because IKewlev. Geller fe Jansenl (120041 ) 
estimated the uncertainty to be ~ 30%, we add this uncer- 
tainly to the observational error as o% ED — a 



[OII]SED 



+ 



{^EWlZ'plte) 2 - 

As discussed in iKawara et al.l (|2010l ). the D(CaII) in- 
dex, which is the ratio of the Call to HI absorption lines in 
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equivalent width, is useful to measure the age of the rela- 
tively old stellar populations. XCaii is defined as: 

XCaii = (D(CaII) obs - D(CaII) temp i at e) 2 /cr% {C aii) (14) 

where (Jo^Caii) is the observational error. 

The luminosity absorbed by dust should mostly be re- 
emitted in the far-IR. If the energy source of the galaxy 
is dominated by star formation, the absorbed luminosity 
should match with the far-IR luminosity. We include this 
effect as follows: 



2 

Xfi 



[Fir -fa&sorfeed) i&ir 



(15) 



where Fi r , Fabsorbed, and a% r are the observed flux in 8- 
1000pm, the absorbed flux in the UV to near-IR, and the 
observed error of Fi r , respectively.lt is noted that Fabsorbed is 
calculated for a given b which is determined by dx%ED/db = 
0. 



6 RESULTS 

First of all, the Calzetti extinction models, which have a 
uniform foreground screen, yield the best fit-model having 
the reduced x-square xi = 5.1 with xfon] = 0-0, XCaii = 
0.3, and Xfir ~ 0-9- This model has a parameter set of the 
optical depth Text (0.3pm) = 5.3, the age t = 2.0 x 10 8 yr 
with the star formation timescale of t s f r = 6.3 x 10 7 yr and 
the metallicity in stellar atmosphere Z — IZq, EW([OII]) 
= 29 A, D(CaII) = 0.09, and F absorbed = 6.6 x 10" 13 ergs 
cm _2 s _1 . This is our reference to which the goodness of our 
other models will be judged. 

The left column panels of Figure|3]plot xt as a function 
of Ti n l{ji n + t sc ) for the uniform foreground screen mixed 
with internal dust. Tj n and t sc are T ex t for internal dust and 
screen dust, respectively. For the pure internal dust geom- 
etry, n n /(nn + t sc ) = 1, while the pure foreground screen 
geometry, Tinjijin + t sc ) = 0. The scattering properties are 
assumed to be the no-scattering, the isotropic scattering (g 
= 0) , or the forward-only scattering (g = 1). In case of no 
scattering, forward scattering is assumed for internal dust. 

The pseudo MW extinction curve, which is made by re- 
moving the 2175 A feature from the MW extinction curve, 
has a smaller xt than the MW extinction curve, suggest- 
ing the 2175 A feature is very weak or unrecognizable in 
SST J1604+4304. Because of the large xl valued relative 
to the Calzetti and SN extinction curves, the MW extinc- 
tion curve is unlikely regardless of the presence of the 2 175 A 
feature. The pure internal dust geometry can be ruled out 
regardless of the choice of extinction curves. Another diffi- 
culty of the pure internal geometry is that too much dust 
is required to reproduce the observed SED (i.e., T ext = 70 - 
200). The estimated dust mass is far greater than that ob- 
served by a factor of 10 or even more. The general trend 
is that xt becomes smaller as Ti n /(Ti n + t sc ) approaching 
to zero, indicating that the extinction due to internal dust 
is a small fraction of the total ext inction as already foun d 
in highly reddened HII regions fi.e. lNatta fc Panagia|[l983 ). 



5 The approximate values of the increment of xt are 1.14, 1.75, 
and 2.36 for 68%(l(x), 95%(2(x), and 99.6%( 3ff) confidence where 
the degree of freedom is v = 12 | Avnii 1 9761) . 
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Figure 5. Sample illustration of the total optical depth Ti„ + r sc 
at 0.3|xm (upper left), the metallicity Z in units of the solar value 
(bottom left), the time after the onset of the first star formation 
t(yr) (upper right), and the e-folding timescale of the star forma- 
tion t s f r as a function of the ratio of the optical depth of internal 
dust to internal+screen dust Ti n /(Ti n + t sc ). The isotropic scat- 
tering (g = 0)is assumed. Solid lines represent SN dust with 8 - 
30 M , dash SN dust with 8 - 40 M Q , dash dot MW dust, and 
dash dot dot pseudo MW dust without the 2175 A feature. 



The two SN extinction curves have Xv l ess than that of the 
Calzetti extinction curve (i.e., xt = 5.1) for a typical range 

Of T in I (Tin + T ac ) = 0.0 - 0.7. 

The middle and right column panels of Figure [4] plot xl 
as a function of the number of clumps along the line of sight 
for the clumpy foreground screen; Pure clumpy foreground 
screens (i.e.,Ti n /(Ti n + t sc ) =0) are shown in the middle col- 
umn panels and foreground screens with internal dust having 
Tin/(Tin + t sc ) = 0.4 in the right column panels. The MW 
extinction curves with and without the 2175 A feature have 
larger xt than that of the Calzetti extinction curve except 
for clumpy screens with the average number of clumps of iV 
= 4 (see the upper right panel). The overall trends are that 
small xt values (i.e. xl < 5.1) are found in larger N. Uni- 
form dust screens have been explicitly or implicitly assumed 
in many cases in the literature. These figures indicate that 
uniform dust screens are good approximations to the real 
distribution of dust. The two SN extinction curves have xt 
smaller than that of the Calzetti extinction curve for large 
N. 

Figure [5] gives typical features of parameters as a func- 
tion of Tin/ {Tin + t sc ) , taking isotropic scattering models as 
an example. Those parameters are the total optical depth 
Ttotai (0.3pm) (= Tin + t sc at 0.3pm), the metallicity Z, the 
time after the onset of the first star formation t, and the 
e-folding timescale t s f r - Because the effect of scattering re- 
duces the effective optical depth r e ff as seen in equations 
(7) to (9), scattering dust screens require a larger opti- 
cal depth than no-scattering screens; taking models with- 
out internal dust , T ex t (0.3pm) ~ 5 for the no-scattering 
screen, T ex t (0.3pm) ~ 8 for the isotropic scattering, and 
Tea:* (0.3pm) ~ 15 for the forward-only scattering. The op- 
tical depth of the SN 8-40 Mq extinction curve requires 





larger optical depths than those of the other curves. This is 
because albedo is greater in the SN 8-40 Mq curve than the 
others (see Table [TJ. Z is smaller on the SN curves than the 
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Figure 4. The left column panels show xl( u = 12) as a function of Ti n /(Tj n + t sc ) for no scattering (top), forward scattering (middle), 
and isotropic scattering (bottom), where T; n is the optical depth for internal dust and r sc for screen dust. The uniform screen is assumed. 
In case of no scattering, forward scattering is assumed for internal dust. Solid lines represent SN dust with 8-30 Mq, dash SN dust 
with 8-40 Mq, dash dot MW dust, and dash dot dot pseudo MW dust without the 2175 A feature. The middle column panels are 
same as the left panels, but as a function of the number of clumps in the line of sight contained in the screen for no internal dust (i.e., 
Tin I (^in +t sc ) = 0.0). The right column panels are same as the left panels, but as a function of the number of clumps in the line of sight 
contained in the screen for Ti n /(ri n + r sc ) = 0.4. xl = 5-1 results if the Calzetti extinction curve is applied. 



MW curve. On the other hand, t is older on the SN curves 
than the MW curves. The older t on the SN curves can be 
attributed to the fact that the slope of the SN curves are 
flatter than that of the MW curve. To reproduce the red- 
ness of the SED, the steep MW curve requires bluer and 
younger stellar populations. 

Figure shows an example of model SEDs along with 
the observed SED. It is noted again that the longest wave- 
length data (i.e., IRAC 8.0u.m) is not used for x3-fitting. 



The sample model is obtained by assuming the SN 8-40 Mq 
extinction curve with the isotropic scattering without inter- 
nal dust. The fitted parameters are T sc (0.3urn) = 10.5, t — 
3.2 xlO 8 yr, t sfr = 1.3 xlO 8 yr, and Z = 1Z Q . This model 
has xl = 3.7 with xfoii] = °- 4 > Xcaii = 0.6, and x% r = 2 A 
providing a smaller xl than those with the Calzetti extinc- 
tion curve. Comparing with xl = 5.1, xfoii] = Xcaii = 
0.3, and Xfir ~ 0-9 for the best-fit model with the Calzetti 
extinction curve, the reduction in xl m this SN model is 
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Figure 6. Sample illustration of model SEDs. The illustrated 
model is obtained by assuming the SN 8-40 Mq extinction curve 
with the isotropic scattering without internal dust. This model 
has xt = 3-7, Tsc = 10.2 at 0.3|xm, t = 5.1 xlO 8 yr, t sfr = 2.5 
XlO 8 yr, and Z = IZq. The top panel shows the model spectrum 
along with the observed SED. The upper curve is the unobscured 
model spectrum. The bottom panel shows xi as a function of t 
and r sc at 0.3u.m. 



attributed to the reduction in Xsed> indicating that the 
reddening by the SN extinction curve reproduces the SED 
better than the Calzetti extinction curve. 



7 SN II EXTINCTION CURVE 

As shown so fa r ; the SN extinction curves synthesized by 
iHirashita et al. (2005) reasonably reproduce the observed 
characteristics of ULIRG SST J1604+4303. The goodness of 
fit is comparable to or better than those with the Calzetti 
extinction curve. The SN extinction curves must be useful 
to study young dusty galaxies not only at high-redshift but 
also at modest-redshift (z ~ 1). We thus present the details 
of the SN extinction curves from the UV to submillimeter. 

Figure shows « a i>s(A) of the SN dust along with that 
of the MW dust, where K a 6s(A) is the mass absorption cross 
section in units of cm 2 g _1 . Table [1] gives K, ext (X) the mass 
extinction cross section, ft a f, s (A) the mass absorption cross 
section, and albedo for the SN dust; SN 8-30 Mq dust in 2-4 
columns and SN 8-40 Mq dust in 5-7 columns . It is noted 
that the o ptical constant of cry stalline SiC>2 is used in the 
models bv lHirashita et al l (2005). As a result, the SN curves 
have spiky features between 9 and 21 (j.m. The contribution 



of crystalline Si02 to the SN extinction curve is a few per- 
cent at most in the mid- to far-IR range except the spikes at 
8.7, 12.25, and 20.5 [im. We smooth the calculated mass ex- 
tinction coefficient of Si grains around 1 fim to remove small 
and somehow artificial ripple features arising from large Si 
grains. 

As shown in Figure [3] the SN extinction curves are flat- 
ter than the others: R v (= A v /E B -v) is 7.8 for the 8-30 
M P) curve and 5.8 for the 8-40 M m , while R v = 4.05 for the 



Calzetti curve and Rv = 3.1 for the MW curve. In other 
words, for a given amount of reddening, the extinction is 
greater in the SN curves than the others. As suggested in 
Figure [7] the silicate feature around 9.7 urn is much less in 
the SN dust than in the MW dust. This feature could be 
used to distinguish the SN extinction curve from the others 
by means of mid-IR observations in the future. 

Another striking feature of the SN curve is the behav- 
ior of /tats (A) from 40 to 1000 urn. Assuming a simple an- 
alytic form of « a 5s(A) oc v 13 , we have J3 ~ 1.6 for the SN 
dust, while /3 ~ 2 for the MW dust. Single-temperature 
dust spectral distributions of the form B v (Tdust), are 
widely used to fit observed SED in the far-IR and sub- 
millimeter, where B v {Td U st) is the Planck function. How- 
ever, fitting only a few datapoints results in a significant 
correl ation between and Tdust that c an account for the 
data (|Blain. Barnard, fc Chapmanll2003T ): higher Tdust with 
smaller /3; Tdust derived with the SN dust is 13% higher than 
that with the MW dust. Assuming the SN dust, we have 
Tdust = 38.3 ± 1AK for SST J1604+4304 which is a com- 
mon characte ristics o f 104 I RAS bright galaxies observed 
at 850 [im bv lDunnd <|2000h . where the best fit was found 
for T du st = 35.6 ± 4.9 K and (3 = 1.3 ± 0.2. 



8 DUST IN SST J1604+4304 

Using Tdust = 38.3 ± 1AK combined with the observed 
160(im flux density and K a f, s (80um) = 22cm 2 g _1 in Ta- 
ble Q] implies that the dust mass M du st in SST J1604+4303 
is 2.2 ± 1.1 x 10 s Mq. Consider a spherical shell with a ra- 
dius of r and a thickness of 5r, and suppose that internal 
dust fills inside the sphere with r and screen dust in the 
shell. The mass of the internal dust can be simply approx- 
imated to Min = 4iv(Ti n / K, e xt)r 2 /3 and that of the screen 
dust to Msc = 4ir(Tsc/K. ext )r 2 , thus M in /M sc = (T in /r sc )/3. 
The best model for the isotropic scattering uniform screen 
without internal dust has r el ((0.3nm) = 10.5. In this case, 
the inferred radius of the spherical shell is about 3.3 kpc or 
0.4 arcsec which is comparable to the optical size of SST 
J1604+4303. 

If the star formation history of SST J1604+4303 is 
described by t = 3.2 xlO 8 yr, t t fr = 1.3 xlO 8 yr, and 
Z — IZq, then the luminosity to mass ratio is 7.6, im- 
plying a stellar mass of 2.4 xlO 11 Mq. During the last 3.2 
x 10 8 yr, ~ 2.0 x 10 9 SNe II exploded in this galaxy. If the 
dust was mainly supplied by SN explosions, then we have a 
dust production rate of ~ 0.1 Mq/SN by simply dividing 
the dust mass by the number of SNe II. 

Although xt is reduced by adopting the synthesized SN 
extinction curves, the goodness of fit is far from the perfect 
fitting(i.e.,x 2 ~ 1). This may come from various reasons. 
For example, the models for dust formation in SNe II are 
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Figure 7. Plot of re a f, s (A) the absorption cross section per dust 
mass in units of cm 2 /g. Note that spiky features between 9 and 
21 urn on the SN curve are mostly attributed to crystal Si02 and 
MgO; spiky features of crystal SiC>2 occur at 8.7, 12.25, & 20.5 
um while a spike MgO feature at 16.5 \J.m 



still coarse and the contribution from the dust formed in 
AGB stars needs to be included, or the stellar population 
models are too simple and the real stellar population may 
be more complex and consist of stars with different ages with 
diff erent metallicit i es. 

iNozawa et al.l l|2007l ) t ake into account the eff ect of re- 
verse shock destruction and lHirashita et al.l (|2008l) calculate 
the extinction curve based on their result. The extinction 
curves after r everse shock de s tructi on are generally flatter 
than those of iHirashita et all (|2005l ). since small grains are 
selectively destroye d. In this paper, we sim ply adopt the ex- 
tinction curves by IHirashita et all (|2005h . In other words, 
the effect of reverse shock destruction is not taken into ac- 
count. Indeed, the flat extinction curves after reverse shock 
destruction are not consistent with the large reddening re- 
quired for SST J1604+4304. This implies that the reverse 
shock destruction may not be efficient or that the small grain 
producti on by, for example, s hattering in warm ISM is oc- 
curring l|Hirashita et al.| [2010). in which case the extinction 
curve after shattering is similar t o that used in this pap er. 
Moreover, the extinction curve by IHirashita et all (|2005l ) is 
supported by observations in the sense th at it is consis- 
tent w ith the extinction curve derived by iMaiolino et al.l 
p004) for a h i gh-z B AL QSO at z = 6.2. Note that the 
Maiolino et all l|2004l) 's extinction curve is also consistent 
with a high-z (z = 6.3) gamma-ray burst (|Stratta et al.l 
120071) . 

We have shown that a flatter extinction curve (or 
a larger R v ) than the MW curve fits the SED of SST 
J1604+4304 better. A flat extinction curve may also be re - 
alized if dust growth occurs in dense clouds (|Drainei r2003). 
Thus, it is worth investigating whether or not a flat extinc- 
tion curve predicted from the grain growth fits the observed 
SED of SST J1604+4304. This is left for future work. 



9 CONCLUSIONS 

We have applied the SN extinction curves to reproduce the 
observed properties of SST J 1604+4304 which is a young IR 



galaxy at z ~ 1. The SN extinction curves used in this work 
was obtained from models of unmixed ejecta of SNe II for 
the Salpeter IMF with a mass range from 8 to 30 Mq or 
8 to 40 Mq. These curves are compared with the MW and 
Calzetti extinction curves. The data used in the \ 2 fitting 
methods are the fluxes in 13 photometric bands from B to 
IRAC 5.8um, the equivalent width of the [OII]A 3727, the 
D(CaII) index, and the IR flux in 8-1000um. 

The effect of dust distributions on the attenuation 
of starlight was investigated by performing the x 2 fitting 
method with various dust distributions. These are the com- 
monly used uniform dust screen, the clumpy dust screen, and 
the internal dust geometry. In addition to these geometries, 
we have investigated three scattering properties, namely, no- 
scattering, isotropic scattering, and forward-only scattering. 
Judging from the \ 2 values, the uniform screen models with 
any scattering property provide good approximations to the 
real dust geometry. Internal dust is inefficient to attenuate 
starlight and thus cannot be the dominant source of the ex- 
tinction. 

We have found that the SN extinction curves reproduce 
the data of SST J1604+4304 comparable or better than the 
Calzetti extinction curve. The MW extinction curve is too 
steep and is not satisfactory agreement with the data un- 
less several dusty clumps are in the line of sight. This trend 
may be explained by the abundance of SN-origin dust in 
these galaxies; SN dust is the most abundant in the young 
IR galaxy at z ~ 1, abundant in local starbursts, and less 
abundant in the Galaxy. If dust in SST J1604+4304 is domi- 
nated by SN dust, the dust production rate is ~ 0.1 Mq/SN. 

We provide the extinction cross sections of the SN dust 
along with albedo from UV to submillimeter wavelengths in 
Table Q] and on-line material table. 
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Table 1. Extinction cross sections of the SN II dust models 



A(|xm) 




bN 8-30 Mq 






bN 8-40 




/ 2 — 1\ 
K ext (cm z g 1 ) 


/ 2 — 1\ 


albedo 


/ 2 — 1\ 


/ 2 — 1\ 


albedo 


6.67c-02 


4.69e+04 


2.34c+04 


5.01e-01 


4.86e+04 


2.37c+04 


5.11c-01 


9.12e-02 


4.73c+04 


2.34c+04 


5.05e-01 


4.95e+04 


2.44e+04 


5.07c-01 


1.00e-01 


4.71c+04 


2.30e+04 


5.12e-01 


4.95e+04 


2.44c+04 


5.07c-01 


1.20e-01 


5.07c+04 


2.57c+04 


4.93e-01 


5.23e+04 


2.64c+04 


4.96c-01 


1.40e-01 


5.09e+04 


2.43c+04 


5.23e-01 


5.25e+04 


2.47c+04 


5.29e-01 


1.60e-01 


4.68c+04 


2.03e+04 


5.67c-01 


4.97c+04 


2.12e+04 


5.74c-01 


1.80e-01 


4.04c+04 


1.67c+04 


5.87c-01 


4.54c+04 


1.61c+04 


6.46c-01 


1.90e-01 


3.87c+04 


1.57c+04 


5.94e-01 


4.29c+04 


1.40c+04 


6.74c-01 


2.00c-01 


3.79c+04 


1.50c+04 


6.04e-01 


4.14e+04 


1.23c+04 


7.02e-01 


2.10c-01 


3.77c+04 


1.46c+04 


6.12e-01 


4.04e+04 


1.13e+04 


7.19c-01 


2.20e-01 


3.77c+04 


1.45e+04 


6.15e-01 


3.96e+04 


1.08c+04 


7.26c-01 


2.30e-01 


3.77c+04 


1.46c+04 


6.14e-01 


3.88e+04 


1.07c+04 


7.24e-01 


2.50e-01 


3.75c+04 


1.44c+04 


6.16e-01 


3.72e+04 


1.04e+04 


7.21e-01 


3.00c-01 


3.40c+04 


1.39e+04 


5.92e-01 


3.10c+04 


9.80c+03 


6.84e-01 


3.55c-01 


3.17c+04 


1.32e+04 


5.82e-01 


2.73e+04 


9.00c+03 


6.70c-01 


3.64e-01 


3.14c+04 


1.31e+04 


5.81c-01 


2.68e+04 


8.89e+03 


6.68c-01 


4.41e-01 


2.80c+04 


1.20c+04 


5.71c-01 


2.26c+04 


8.15e+03 


6.40c-01 


4.69c-01 


2.69c+04 


1.16c+04 


5.69c-01 


2.14c+04 


7.84c+03 


6.34e-01 


5.47c-01 


2.48c+04 


1.08c+04 


5.65c-01 


1.93c+04 


7.30c+03 


6.22c-01 


6.17c-01 


2.21c+04 


9.88e+03 


5.53e-01 


1.70c+04 


6.64e+03 


6.10c-01 


6.42e-01 


2.10c+04 


9.66e+03 


5.40c-01 


1.61e+04 


6.51c+03 


5.96c-01 


6.49c-01 


2.08c+04 


9.60e+03 


5.39e-01 


1.60c+04 


6.47e+03 


5.95e-01 


7.00c-01 


2.03c+04 


9.21c+03 


5.46e-01 


1.55c+04 


6.19c+03 


6.00c-01 


7.48c-01 


1.95c+04 


8.90e+03 


5.45e-01 


1.49e+04 


5.98c+03 


5.99c-01 


8.02c-01 


1.86c+04 


8.55c+03 


5.40e-01 


1.41c+04 


5.73e+03 


5.94e-01 


8.66e-01 


1.75c+04 


8.16c+03 


5.34e-01 


1.32c+04 


5.46e+03 


5.86c-01 


8.93e-01 


1.70c+04 
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APPENDIX A: TWO STREAM 
APPROXIMATION FOR RADIATION 
TRANSFER IN A DUSTY SLAB 



The two stream approximation in radiation transfer was ap- 
plied to a sph e rically symmetric dust envelope surrounding a 
star bv lCodel 1 19731 ). Here we consider a plane-parallel dust 
slab. 

We measure distance z outward for a plane-parallel slab, 

then 



- ^ = j v - K v I v (z, At) + t^- / p(fi, n')I v (z, [i)dw' (Al) 

P dz 47r J An 

where I v (z,y?j is the specific intensity, /i = cos,9 the direc- 
tional cosine, v the frequency, k„ the mass extinction coeffi- 
cient, cr„ the mass scattering coefficient, p the mass density 
of dust, j v the ratio of the stellar emission per unit dust 
mass. The phase function fj,') is defined by 



p(n, n')duj' = 4tt or - / p(fi, fi')d(i' = 1 (A2) 

4tt 1 J-l 

By introducing the optical depth r v which is defined by 
dr u — —pK, u dz, the radiation transfer equation is reduced to 



dr u 



1» 



p(\i,\j!)I v (t, n')dij! (A3) 



where 7„ = ov/fi^ is the albedo. In what follows, we omit 
the subscript v for simplicity. 



Al Slab illuminated at the back side 

Consider a slab which contains no emitters and illumi- 
nated at the back side, i.e., j„ = 0. In the two stream 
approximation, the specific intensity n) is represented 
by the specific intensities at two specified directions at ±/x; 
7 + (r) = I(t,h) and J"(r) = I(r,-fi) for < fj, < 1. 7 + (r) 
and I~(t) are assumed to be independent of n over the up- 
ward and the dow nward hemis pheres, respectively. 

According to I Code! l|l973h . we introduce the forward 
and backward scattering coefficients a + and a~ which are 
defined by: 

— = 7;[ v{n,n')dii and — = \ [ p(/j,, (A4) 
a 2 Jo a 2 Jo 
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respectively. Note that a + + a = a. A modified asymmetry 
factor g is defined by 



1.00 



g=- (A5) 

a 

note that g = 1 for forward scattering, g — — 1 for backward 
scattering, and g = for isotropic scattering. Then using 
the relations a + /a = (1 + g)/2 and a~ /a = (1 — g)/2, the 
equations (A3) is reduced to: 



W 



0.10 



dl+(l 



I + (r) 



f (l±l)/+( T )- 7 



(1-3) 



J-(t) (A6) 



- = /-(r) - 7^/"(r) - 7^/ + M (AT) 

Here we replace the optical depth along the z-axis with 
the optical depth along the line of sight dText — dr/fi and 
introduce new variables x = I + — I~ and y = I + + I~ , the 
above equations become simply: 



dx , . dy , . 
= (1-7) y and = {l-jg)x 



(A8) 



diT ex t dText 

Then the solutions of I (r ex t) and I~(r ex t) are given 



by 



1 



I + {r ext ) = - P+ l)Aexp (area*) - (/3 - 1)B exp (-ttT elt )] (A9) 



/" (Tarf) = - [(/3 - 1)A exp (ar ell ) - (/? + 1)B exp (-aText)KAlO) 

where a 2 = (1 - 7 )(1 - 73) and /? = [(1 - 7ff)/(l - 7)] 1/2 . 
Using the boundary conditions of 7~(0) = and I + (Ti ext ) = 
Io for determing the constants A and B, where Ti e xt is the 
optical depth of the slab along the line of sight, then the 
emergent intensity I + (0) is given by 



7+(0) 



2/3/q 



(j3 + 1) exp (ariext) + (/?-!) exp (-aTiext) 



(All) 



A2 Slab containing emitters 

We now consider a slab which contains emitters ,i.e., stars. 
The slab is illuminated by these internal stars only. Then 
the radiation transfer equations for I are simply given by 



dl+(i 

1 dT 



= J+(r) _ 7 (l+s) 7 + (r) 
-7^J-(r)-A 



dl-{r) 

dT 



= /- (r )_ 7 C!+£lj-( T ) 



(A12) 



(A13) 



where A = ] u /k v . Then x = I + — I andy = 7 + + 7 can be 
determined by solving the following simultaneous equations; 



0.01 



0.01 0.10 1.00 10.00 100.00 



''ext 



Figure Al. Plot of the apparent extinction e; n as a function of 
the optical depth along the line of sight Text- The geometry is 
the internal dust with 7 = 0.65 and g = (spherical scattering). 
The solid line represents the solution of equation (A19) for A 
= const against Text- The other lines represent the absorption 
approximated by e; n = (1 — e~ T "ff ) I Teff ; the dotted line for 
T e// = (! ~ T)Text,th.e dashed line for T e // = (1 — 7 )°- 5 r e:Et , 



the dash dot line for t. 



// 



(l- 7 )°- 75 T e . t . 



The intensities I (Text) are given by 

I + {T ext ) = \(P + 1) [A - So (Test)] exp (oTest) 
-\{P-1)[B- El (T ex t)] eX P (-QT e;ct ) 

I' (Text) = \{fi— 1) [A- 7J (T ext )] exp (otText) 
-\(P+1)[B- El (T e:c t)] CXp (-QTeri) 

where 

Eo(T e xt) 

7^1 (Text) 



Jo'** H T Lxt) exp (-OTeat) ^T^t and 

= /J""*' A(Text) exp (aText) dText 



(A15) 



(A16) 



(AIT) 



Given the boundary conditions of 7 (0) = I + (Ti ex t) = 
0, the emergent intensity 7 + (0) is given by: 

T+/n\ - oa O 3 + l)-^0 (next) exp (QTlext) - (J3 - 1)^1 (next) exp (-OTlggt) , . 1 ^ 

W — , a ; -.no 7 \ 77, 77o 7 \ 

(j8 + l) 2 exp (aTiext) — (p - 1) exp (-QTiext) 

If A = const against Text, 7 + (0) — > Ati s // for Tiext — ► 0. 
Thus, the apparant extinction e^ n as defined in equation (1) 
is: 

2/3 03 + l)[exp (c*T le xt) - 1] - (/? - 1)[1 - exp (-an e xt)] 



aTl ex t 



(13 + l) 2 exp (aTiext) - (P - !) 2 exp (-aT lex t) 



<A19) 



Figure IA1I compares equation (A19) with the approxi- 
mations by e ln = (1 - e~ T "ff)/Tef f for Teff = (1 - 7) r erf, 
fe// = (1 ~ 7)°' 5 ^t, and r e// = (1 - 7 ) a75 r e;c t. The ge- 
ometry is internal dust with 7 = 0.65 and g = 0. As shown 
in Figure [ATI Teff = (1 — j) ' 75 T e xt approximates to equa- 
tion (A19) better than T e f / = (1 — 7 )°' 5 T e xt for the isotropic 
scattering. 



dx 

dText 



= (1 — 7) y — 2A and 



dText 



= (l-7S)z. (A14) 



